Virus shapes and buckling transitions in spherical shells 
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We show that the icosahedral packings of protein capsomeres proposed by Caspar and Klug for 
spherical viruses become unstable to faceting for sufficiently large virus size, in analogy with the 
buckling instability of disclinations in two-dimensional crystals. Our model, based on the nonlinear 
physics of thin elastic shells, produces excellent one parameter fits in real space to the full three- 
dimensional shape of large spherical viruses. The faceted shape depends only on the dimensionless 
Foppl-von Karman number 7 = YR 2 /k, where Y is the two-dimensional Young's modulus of the 
protein shell, k is its bending rigidity and R is the mean virus radius. The shape can be parame- 
terized more quantitatively in terms of a spherical harmonic expansion. We also investigate elastic 
shell theory for extremely large 7, fO 3 < 7 < fO 8 , and find results applicable to icosahedral shapes 
of large vesicles studied with freeze fracture and electron microscopy. 
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I. INTRODUCTION 

Understanding virus structures is a rich and challeng- 
ing problem with a wealth of new information now 
becoming available. Although traditional X-ray crystal- 
lography still allows the most detailed analysis Q , three- 
dimensional reconstructions of icosahedral viruses from 
cryo-electron micrographs are now becoming routine Q . 
Electron microscope images of many identical viruses in 
a variety of orientations are used to reconstruct a three- 
dimensional image on a computer, similar to CT (com- 
puted tomography) scans in medical imaging. There are 
now, in addition, beautiful single molecule experiments 
which measure the work needed to load a virus (bacte- 
riophage 4>29) with its DNA package 0. The aim of 
this paper is to explore the elastic parameters and phys- 
ical ideas which determine the shapes of viruses with an 
icosahedral symmetry, using the theory of thin elastic 
shells 

The analysis of approximately spherical viruses dates 
back to pioneering work by Crick and Watson in 1956 0, 
who argued that the small size of the viral genome re- 
quires identical structural units packed together with an 
icosahedral symmetry. These principles were put on a 
firm basis by Caspar and Klug in 1962 0, who showed 
how the proteins in a viral shell (the "capsid" ) could be 
viewed as icosadeltahedral triangulations of the sphere by 
a set of pentavalent and hexavalent morphological units 
("capsomers" ). The viral shells (there can also be an 
outer envelope composed of additional proteins and mem- 
brane elements from the host cell) are characterized by 
a pair of integers (h, k) such that the number of mor- 
phological units is N = 10(h 2 + hk + k 2 ) + 2. To get 
from one pentavalent capsomer to another, one moves 
h capsomers along a row of near-neighbor bonds on the 




FIG. 1: A right-handed (3,f) triangulated net (icosadeltahe- 
dron) used to describe virus structure. The (f ,3) structure is 
left-handed. 



sphere, turns 120 degrees and moves another k steps. 
Euler's theorem relating the number of vertices, edges 
and faces of a spherical triangulation insures that the 
number of capsomers in five-fold environments is exactly 
12 A simple icosahedron of 12 morphological units 
corresponds to (1,0) while soccer balls and Cgo fullerene 
molecules are (1,1) structures with N = 32 polygons. 
A (3,1) icosadeltahedron N = 132 is shown in Fig. ^ 
The polyoma virus (SV40) is a (2,1) structure with 72 
capsomeres, while the much larger adenovirus and her- 
pes simplex virus are (5,0) and (4,0) structures with 252 
and 162 morphological units, respectively. Structures like 
that in Fig. ^ with h and k nonzero and h ^ k are chiral. 

Note that the relatively small polyoma virus (diameter 
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FIG. 2: The polyoma virus (a) is approximately spherical, 
while the larger adenovirus (b) is more faceted (not to scale). 
Images from Ref. 

440 Angstroms) is round (see Fig. |2Ji), while the much 
larger herpes simplex virus (diameter 1450 Angstroms) 
has a more angular or faceted shape |S( (see Fig. |2t>). 
Faceting of large viruses is in fact a common phe- 
nomenon; the protein subunits of different viruses, more- 
over, are very similar (see below) . If these protein assem- 
blies are characterized by elastic constants and a bending 
rigidity p|, we can ask how deviations from a spherical 
shape develop with increasing virus radius, which scales 
roughly as the square root of the number of morphologi- 
cal units. 

In support of the idea that viruses with different overall 
capsid size are composed of nearly identical monomers, 
we note that most viral coat and capsid proteins have 
about the same size, molecular weight, amino acid com- 
position and, most importantly, the same folded struc- 
ture in three dimensions [To| . It is known, moreover, 
that protein structure determines the mechanical proper- 
ties of proteins Hence, the similarity of the pro- 
tein structure of the coat/capsid proteins suggests similar 
mechanical properties. In addition, the presence of the 
same fold of capsid proteins in unrelated viruses (bac- 
terial phages, plant viruses, insect viruses and animal 
viruses) indicate that the fold and its mechanical prop- 
erties are conserved in evolution and could be essential 
for proper virus assembly. 

In this paper, we argue that the faceting of large 
viruses is caused by a buckling transition associated with 
the 12 isolated points of 5-fold symmetry. These sin- 
gularities can be viewed as disclinations in an otherwise 
6-coordinated medium. It is well known that the large 
strains associated with an isolated disclination in a flat 
disk spanned by a trian gula r lattice leads to buckling into 
a conical shape for (l3lll4j 

YR 2 /k > 154, (1) 

where Y is the two-dimensional Young's modulus, k is the 
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FIG. 3: The fivefold disclination in a triangular lattice. In 
viruses, the points would correspond to capsomeres, while 
in lecithin vesicles they correspond to lipid molecules. The 
highly strained flat space configuration is shown in (a). The 
buckled form which arises for 7 > 154 is shown in (b). 

bending rigidity and R is the disk radius. The energy of a 
single 5-fold disclination with "charge" s — 2ir/6 centered 
in a flat array of proteins of size R is approximately 

E 5 « ^s 2 YR 2 . (2) 

OZ7T 

However, above a critical buckling radius Rb ~ 
^/154k/Y , is there a conical deformation (see Fig. OJ) 
such that the disclination energy now only grows loga- 
rithmically with R 

E 5 a (tt/3)k ]n(R/R b ) + ^-s 2 YR 2 . (3) 

OZ7T 

One might expect a similar phenomenon for 12 discli- 
nations confined to a surface with a spherical topology. 
Indeed, the elastic energy for 12 disclinations on an unde- 
formed sphere of radius R is has a form similar to Eq. J5J), 
namely |l5j | 

E w 0.604(s 2 ri? 2 /47r), (4) 

where the sphere radius R now plays the role of the sys- 
tem size. Although it seems highly likely that these 
12 disclinations can lower their energy by buckling for 
large R, the nonlinear nature of the underlying elastic 
theory leads to complex interactions between the re- 
sulting conical deformations. A boundary layer analysis 
of the von Karman equations for bent plates predicts 
anomalous scaling for the mean curvature in the vicin- 
ity of the ridges connecting conical singularities 0, 0] . 
Interesting scaling behavior also arises in the vicinity of 
the apexes of the cones themselves [l^. Another inter- 
esting physical realization of the buckling problem lies 
in the faceting of lecithin vesicles at temperatures suf- 
ficiently low so that the lipid constituents have crystal- 
lized Ea, i3. 

In this paper, we study the ground states of crystalline 
particle arrays with 12 disclinations in a spherical geome- 
try. We find that there are indeed striking manifestations 
of the buckling transition even in the curved geometry of 
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viral capsids or crystalline vesicles. The nonlinear Foppl- 
von Karman equations for thin shells with elasticity and 
a bending energy are solved using a floating mesh dis- 
cretization developed and studied extensively in the con- 
text of "tethered surface" models of polymerized mem- 
branes (2l| . By taking the nodes of the mesh to coincide 
with the capsomers, even small viruses can be handled in 
this way, although any buckling transition will surely be 
smeared out unless Rb/a is large, where a is the spacing 
between these morphological units. Ideas from contin- 
uum elastic theory will, of course, be most applicable for 
vesicles composed of many lipids and for large viruses — 
Viruses with as many as 1692 morphological units have 
been reported |22l |. 

There may be inherent limitations on the size of vi- 
ral capsids that follow from the elastic properties of thin 
shells. Because larger viruses can accommodate more 
genetic material, large sizes could confer an evolutionary 
advantage. If, however, large viruses buckle away from a 
spherical shape, the resistance of the capsid to mechani- 
cal deformation can degrade. As we shall see, a theory of 
buckled crystalline order on spheres also allows estimates 
of important macroscopic elastic parameters of the cap- 
sid shell from structural data on the shape anisotropy. 
Estimates of quantities such as the bending rigidity and 
Young's modulus of a empty viral shell might allow an 
understanding of deformations due to loading with DNA 
or RNA 0. Although some aspects of virus structure 
may be accounted for by the physics of shell theory, we 
should emphasize that other features could be driven by 
the need for cell recognition, avoidance of immune re- 
sponse, etc. 

A summary of our investigations of buckling tran- 
sitions in a spherical geometry (discussed in detail in 
Sec. II) is shown in Fig. 0] As illustrated in Fig. 0] (See 
also Fig. icosahedral shells do indeed become aspher- 
ical as the "Foppl-von Karman number" 7 = YR 2 /k 
increases from values of order unity to YR 2 /n 3> 1. 
The mean square "asphericity" (deviation from a per- 
fect spherical shape) departs significantly from zero when 
YR 2 /k exceeds 154, the location of the buckling transi- 
tion in flat space 131. Fits of buckled viruses or crys- 
talline vesicles to this universal curve would allow an ex- 
perimental determination of the ratio Y/k. More quanti- 
tative information on the buckled shape can be obtained 
by expanding the radius R{9, 4>) in spherical harmonics, 



R{9,< 



(5) 



and studying the rotationally invariant quadratic invari- 
ants allowed for viruses or vesicles with icosahedral sym- 
metry, namely 
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FIG. 4: Mean square asphericity as a function of YR 2 /n 
for many different icosahedral shells. The inset shows a (6,6) 
structure with Foppl-von Karman number 7 = YR 2 Jk ~ 400. 
The arrow marks the location of the buckling transition in flat 
space. 



rameter set of the form {Q6, Q101 Q12 ■ ■ •} could be con- 
sistent with an icosahedral symmetry, all buckled objects 
describable by the theory of elastic shells in fact lie on 
a universal curve parametrized by the value of YR 2 /k. 
Deviations from this curve would presumably describe 
biological features such as the protrusions of the aden- 
ovirus in Fig. [21 

In Sec. II, we describe our theoretical results for discli- 
nation buckling in the icosadcltahcdral spherical shells 
proposed by Caspar and Klug as models of viruses 0. 
The energy, mean-square aspherity, and spherical har- 
monic content of these shells are determined as a function 
of the Foppl-von Karman number, discussed above, 



7 = YR 2 /k. 



(7) 



with I = 0, 6, 10, 12, 16, 18, • • • jig. Although any pa- 



Most viruses have either Foppl-von Karman number 7 < 
150 (implying a close to spherical shape) or 200 < 7 < 
1500 (noticeably buckled). Higher von Karman numbers 
describing objects with very sharp corners cannot be ob- 
tained for viruses with R < 0.2/im composed of finite 
size proteins. Of course, very high von Karman numbers 
are possible for spherical vesicles with crystalline order 
composed of much smaller lipid molecules 0, |2(| • 

We have studied the scaling of the curvature C at the 
creases formed after the shells buckle for large 7. We 
eventually recover the scaling proposed and studied in 
other geometries by Lobkovsky, Witten and collabora- 
tors but only for very large 7, 7 > 10 7 , appro- 
priate for the buckled icosahedral vesicles described in 
Ref. [23. 

In Sec. Ill, we discuss briefly the relevance of our work 
to icosahedral viruses from the library of those whose 
structures have been determined by diffraction methods 
or by cryro-electron microscopy 0, LJ y> Oj ■ For viruses 
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large enough to buckle, the model seems to account well 
for the deviation from the spherical shape using the sin- 
gle adjustable parameter 7. These fits in turn provide 
information about the ratio Y/re of the Young's modu- 
lus to the bending rigidity. We also discuss the possible 
relevance of spontaneous curvature terms and outward 
pressure induced by packaged DNA and RNA in this Sec- 
tion. See Ref. [24| for a discussion of similar issues for the 
shapes of liquid membranes with a spherical topology. 



II. DEFECTS ON CURVED SURFACES 



Topological defects play a very important role in crys- 
talline matter. A particularly common type of defect, 
the dislocation, is largely responsible for the strength of 
materials, and in two-dimensional systems the unbind- 
ing of dislocations may drive the crystal into a hexatic 
phase 25] . Disclinations on the other hand, are much 
less common in the crystalline phase because of their very 
large energy. In quasi- two-dimensional curved surfaces, 
however, the situation may be quite different. The Gaus- 
sian curvature of the surface will "screen" out the strain 
around the defect and thus lower the energy. Moreover, 
when a crystalline surface is bent to form a closed sur- 
face with spherical topology, defects are necessarily in- 
troduced into the lattice. For a triangular lattice on a 
sphere the number of disclinations has to be at least 12. 
More generally, the difference N5 — Nj between the num- 
ber of fivefold disclinations and the number of sevenfold 
disclinations (assuming defects with coordination num- 
bers other than 5, 6 or 7 are absent) is precisely 6 times 
the Euler characteristic \ of the triangulated surface, i.e., 
N$ — Ny = 6\ = 12(1 — <?), where g is the genus or number 
of handles. Thus for shapes such that g = (the sphere) 
and torii with extra handles (g > 2), the ground state 
must necessarily contain disclinations. (A simple torus 
has genus g = 1, so there is no topological necessity for 
defects.) 

In this section we study closed surfaces of spherical 
topology. The 12 disclinations present in the crystalline 
lattice structure can be expected to dominate the en- 
ergetics and affect the overall shape of the structure. 
The repulsive interaction between the 12 disclinations 
will favor an arrangement which maximizes their separa- 
tion. Absent an instability toward grain boundaries |15| | , 
this leads to a configuration with icosahedral symmetry, 
where the disclinations sit at the vertices. We show below 
that, as a result of competition between strain and bend- 
ing energies, these structures may undergo a buckling-like 
transition (smeared by finite size effects) from a smooth 
round shape to a sharply faceted shape as the size or 
elastic constants are varied. 



A. Disclination Buckling on Spheres 

We assume a thin shell described by a continuum elas- 
tic theory, with energy H = H s + Tib [5j, including both 
an in-plane stretching energy, 



K S = \J dS (2fxu 



Xu lk) 



(8) 



where Uij is the strain tensor, \x and A are the 2D Lame 
coefficients, and a bending energy, 



J dS (kH 2 + 2k g K) , 



(9) 



where re is the bending rigidity, kg the Gaussian rigidity, 
H and K the mean and Gaussian curvatures, respec- 
tively. (If i?i and R2 are the principal radii of curvature, 
H = 1/Ri + 1/R 2 and K = l/R^ Q.) For a closed 
surface with fixed topology the Gaussian curvature inte- 
grates to a constant (provided that kg is constant) and 
will henceforth be dropped, as it will have no influence 
on the shape. Instead of the Lame coefficients we will 
use the 2D Young's modulus Y and the Poisson ratio v, 
which are given by 



Y 



+ X : 



A 



2/i + A 



(10) 



Taking the variation of TL with respect to coordinates pa- 
rameterizing the surface leads to the Foppl-von Karman 
equations which are highly nonlinear and whose solution 
even in simple geometries is very difficult [f|. 

The energy of disclinations on flexible crystalline mem- 
branes was studied by Seung and Nelson in Ref. . As 
discussed in the Introduction, for a thin flat plate of fi- 
nite radius R with an isolated 5-fold disclination at the 
center, the energy grows quadratically with the radius, 
-Eftat ~ AYR 2 , where A ps 7r/288 is a numerical con- 
stant. If the disclination is allowed to buckle out of the 
plane this energy is reduced, and grows logarithmically 
for large R with a coefficient proportional to re. In the 
inextensional limit Y — > 00 the problem simplifies consid- 
erably and the energy is -Ebuckied — Bk In (R/a), where 
B w 7r/3 and a the lattice constant. Thus, for small 
plates the flat solution is lower in energy, whereas for 
large plates the buckled solution wins. An instability sep- 
arating these two regimes occurs when E^t ~ ^buckkdj 
or YR 2 /k sa B/A. A detailed calculation in Ref. fia ] 
found the transition at a critical value of YR 2 / k « 154. 

The spherical shapes of icosahedral symmetry we are 
interested in here can approximately be thought of as 
being composed of 12 disclinations, and should there- 
fore undergo a similar transition from flat to buckled. 
Because the surface of the sphere is already curved, and 
hence breaks the up-down symmetry that was present for 
a thin plate, it is not clear that a sharp instability sur- 
vives in this case. However, even if this is the case, we 
might still expect to see remnants of the transition in the 
form of a sharp if not singular crossover. We construct 



below simple estimates of the energies and transitions 
involved. 

The total energy of the closed shell in the vicinity of 
the transition is approximately 12 times the energy of 
a disclination (with radius approximately equal to the 
radius of the sphere), plus contributions from the back- 
ground curvature of a sphere, given approximately by 
87TK + 47r/tG. In the inextensional limit the short distance 
cutoff in the buckled disclination energy was provided by 
the lattice constant. For finite Y the cutoff is determined 
instead by a balance of strain and bending. We may in 
this case approximate the buckled disclination by a flat 
inner region r < R\, with energy AY R 2 and an outer 
buckled region r > i2& with energy Bn\Yi(R/Rb). Mini- 
mizing the sum with respect to Rb gives YR^/k = B/2A, 
and E ~ kB (| + In (R/R b )) ■ 

Because the disclination energy is independent of the 
two-dimensional Poisson ratio and the Gaussian rigid- 
ity drops out, the solution depends only on a single 
dimensionless parameter, 7 — YR 2 /k, which we term 
the Foppl-von Karman number. Note that, if the 2D 
elastic theory derives from a thin shell of finite thick- 
ness d built from a 3D isotropic elastic medium, we 
have 7 = 12 (l — i/f) (R/d) 2 , where z/3 is the three- 
dimensional Poisson ratio, a result independent of the 
3D Young's modulus Y 3 p| . In summary, we then expect 
the energy of the closed shell with twelve disclinations to 
approximately be 




FIG. 5: Numerically calculated shapes with (h, k) indices 
(2,2), (4,4), (6,6), and (8,8) for fixed k = 0.25, and fixed 
spring constant e — 1. The Foppl-von Karman numbers for 
these shapes are 7 « 45, 176, 393 and 694. 



B. Numerical Results 



E 



6^7/76 + D, 
65(1 Tin ( 7 / 76 )) 



D, 



7 < lb 
7 > lb 



(11) 



where 76 = YR 2 /n. The background curvature gives a 
constant contribution D w 47r (2 + kq/k), leading as well 
to a small shift in the disclination energies and hence in A 
and B. For a perfect sphere, e.g., 12A ps 0.604tt/36 [T^ . 

As the disclinations become more sharply buckled the 
whole structure will become more faceted. This leads 
to the formation of ridges connecting the vertices of 
the icosahedron. The energy of similar ridges has been 
studied recently 0, 0], uncovering some remarkable 
scaling relations. As the ridges become sharper upon 
increasing 7, the energy will increase as -Eridge/^ ~ 

1.24a 7 / 3 (YL 2 /k) 1/G , where a is the angle in radians and 
L is the length of the ridge. In this regime the shape 
is very close to an icosahedron with sharp facets, with 
a w 0.365 and L « 1.23-R, and a total of 30 ridges. In 
the limit of very large 7 = YR 2 /k the energy should 
therefore crossover to 



E/k ~ C7 1/6 + const., 7 -> 00, 



(12) 



where C ss 3.8. 

To check these arguments and to calculate more prop- 
erties of the shells we now present numerical calculations. 



For numerical calculations it is useful to consider dis- 
cretized versions of Eqs. © and © 0: 



2 ^ 



Oi-Tjl -of 



and 



Ub = - ^2 (nj - hjf . 

(U) 



(13) 



(14) 



Here (ij) denote pairs of nearest neighbor vertices (which 
we identify with the centers of the capsomers of a virus), 
with positions r^, and (IJ) pairs of nearest neighbor pla- 
quettes of a triangulated surface, with unit normals nj. 
In the continuum limit this model becomes equivalent to 
Eqs. © and JO} with parameters 01 



K 2 K, 



1 



G 3 



(15) 
(16) 



The relation kg — — 4k/3 was calculated by comparing 
the bending energy of a triangulated cylinder and sphere 
with the corresponding continuum expressions, and dif- 
fers from the one used in Ref. 13]- Closed triangular 
surfaces of icosahedral symmetry are constructed for non- 
negative integers (h, k) according to the geometric prin- 
ciples of Caspar and Klug 0, and the minimum energy 
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y=YR 2 /K 

FIG. 7: Mean squared aspherity (see Eq. 1171 1 as a function 
FIG. 6: Total energy. Dotted and dashed lines are fits to of 7 = YR 2 /k for various "triangulation numbers" T = h + 
Eq. 1111 and 1121 . respectively. The arrow indicates the value hk + k 2 . Data from shapes of three different sizes collapse 
of 7b obtained from the fit. onto a single universal curve, indicating that the continuum 



limit has been reached for these sizes. 



configuration is found numerically using a conjugate gra- 
dient method for different values of k. As discussed in 
the Introduction, the integers (h, k) denote the number of 
steps along the two lattice vectors between two neighbor- 
ing disclinations of the structure. Figure |S] shows some 
examples of the resulting shapes. Shapes of varying size 
with "T-numbers" as large as T = h 2 + hk + k 2 m 1500 
are studied. Below we calculate some properties of the 
shells to characterize the shapes quantitatively. 



1. Energy 

We first plot the total energy as a function of the Foppl- 
von Karman number 7 in Fig. |SJ Similar to results for 
an isolated disclination in a disk with free boundary con- 
ditions |13j, the energy crosses over from a "flat" regime 
dominated by stretching energy to a "buckled" regime 
dominated by bending energies. Fitting the functional 
form Eq. (TTJ gives A m 0.005, B sw 1.30, j b m 130, 
which compares quite well with the estimates above. For 
large 7 > 10 7 a crossover to the form given by Eq. i|12|) 
(indicated by the dashed line in the figure) occurs. 



2. Aspherity 



where Ri is the radial distance of vertex i and (R) is the 
mean radius, 



(R) 



1 N 

-Yr 



(18) 



1=1 



The result, which is displayed in Fig. [7| shows a rather 
sharp, but non-singular crossover from spherical shape 
to faceted at roughly 7 = YR 2 /n w 150. The second 
increase around YR 2 /k w 10 7 coincides roughly with the 
sharpening of the ridges, where the asymptotic scaling in 
Eq. ijT^} sets in. Note that the log-linear plot of Fig. 
extends the range of 7 in Fig. 01 by six orders of magni- 
tude. 



3. Icosahedral spherical harmonics 

We first expand the radial density of points on the 
surface in spherical harmonics, 



(19) 



where the density R(9, 4>) is defined by 

R(0,<j>) = ^2R j 5((t>-(l> j )5(cosd-co86 j ), (20) 



As a measure of the deviation from a perfectly spher- 
ical shape centered on the origin we calculate the mean 
squared aspherity, defined by 



(AR 2 ) 
(R) 2 



1 N 

-Y 



(R) 2 



(17) 



and (Rj , 9j , <f>j ) represents the polar coordinates of ver- 
tex j. From the coefficients Q; m , we form the rotation- 
invariant combinations 



Qi - 21 + 1 Y, \Qlm\ 



(21) 
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FIG. 8: Combination of spherical harmonics (see text), which 
fall on a universal curve for large enough shells. Points for 
triangulations denned by (h,0),h — 5, 6, . . . , 25 were used 
to construct the solid curve. For h < 5 deviations due 
to discreteness become large and those points have been 
omitted for clarity. Points labeled 1-6 correspond to 7 ~ 
0.5, 30, 1000, 15000, 8 x 10 6 and 2.5 x 10 8 , respectively. 



For a shape of icosahedral symmetry the Qis are non- 
zero only for I = 0, 6, 10, . . . We plot, in Fig. El 
Qw/Qo vs Qq/Qo, which for large shells should fall on 
a universal curve parameterized by YR 2 /k. Note that 
any point in the (Qe/Qo, <3io/Qo)-pla n e would be con- 
sistent with an icosahedral symmetry. Continuum elas- 
tic theory, however, predicts a universal set of functions 
Qi (7) parametrized only by the Foppl-von Karman num- 
ber. The buckling transition occurs between the points 
labeled 2-4 in Fig. |S| while the crossover to the ridge 
scaling for very large 7 happens between points 4-6. 



4- Curvature 

The curvature C across the midpoints of the ridges con- 
necting the vertices of the icosahedra is plotted in Fig. [5] 
As 7 = YR 2 /k increases through the transition the ridges 
get sharper and the shape becomes more faceted. A per- 
fect sphere would have C = 1/R. However, as seen in 
Fig. El CR saturates to a slightly smaller value, ~ 0.7, 
for very small 7, implying that the shape is not perfectly 
spherical below the buckling transition. In fact, there is 
a weak tendency toward a dodecahedral shape (which is 
the dual to the icosahedron) . The effect is hardly visible 
(cf. Fig. EJ), however. The data for large shells is well 
described by a scaling form 



C = R- 1 F{YR 2 /k) 

depending on the single parameter 7 
arguments 7 — > 00, we find F(^y) — ► 7 



(22) 

In the limit of large 
1 / 6 consistent with 
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FIG. 9: Curvature at the midpoints of the ridges. The points 
labeled 1-6 are for the same values of 7 as in Fig. El The 



dashed line shows the asymptotic scaling behavior ~ 7 



1/6 



the scaling arguments of Lobkovsky et al. |16lll7| . Note, 
however, that 7's well in excess of 10 6 are required before 
one begins to see this asymptotic result of "ridge scaling" . 



III. DISCUSSION 

We have analyzed a model, based on the (highly non- 
linear) physics of thin elastic shells, which may be suit- 
able for describing the shapes of large viruses and of large 
vesicles with crystalline order in the lipids. Application 
of our results to vesicles 0, |2(j seems straightforward, 
once sufficiently precise freeze fracture or confocal mi- 
croscope images become available. Figure ITU1 illustrates 
two highly faceted shapes we found for the large Foppl- 
von Karman numbers or "vK's" which might be relevant 
to the experiments of Ref. [2(j. One complication ne- 
glected here concerns possible phase separation of the bi- 
nary lipid mixtures studied by Dubois et al. in the vicin- 
ity of the 12 disclinations [20|. It would be interesting 
to investigate this effect, although a modest enrichment 
of one lipid species near a disclination could be incorpo- 
rated into a renormalized core energy. As discussed in 
Sec. Ill Bl vK's in excess of 10 6 are required to see clearly 
the interesting scaling predictions of Refs. ^(| and |17| . 
In the remainder of this section, we comment on the rel- 
evance of our work to spherical viruses 0, 0, 0, , where 
the vK's are of order a few thousand or less. 

Ref. Q compiles cryro-electron micrographs and other 
data on approximately thirty different viruses, arranged 
in order of increasing size. These images highlight the 
trend that small viruses are round and larger viruses 
are more faceted. We view faceted viruses as the re- 
sult of 12 simultaneous buckling transitions, centered on 





FIG. 10: Two shapes for the large Foppl-von Karman num- 
bers 7 = 15600 and 8 x 10 6 (points 4 and 5 in Fig. 
illustrating the sharpening of the ridges in the ridge scaling 
regime. 



12 disclinations, similar to the buckling of an isolated 
disclination centered on a disk with open boundary con- 
ditions 01 . The spherical packing of the protein cap- 
somers in viruses not only forces in twelve disclinations 
(which we assume reside at the vertices of an icosahe- 
dron) 7], but also breaks the up/down symmetry of a 
disk with respect to the direction of buckling. Hence, 
we expect (and find numerically) that the sharp buck- 
ling transition with increasing size in Ref. is smeared 
out. As shown in Fig. ^] for bacteriophage HK97 (with 
72 capsomers), good one parameter fits in real space to 
the full three dimensional shape of spherical viruses are 
possible. Our best fit for Foppl-von Karman number of 
this mature form of HK97 is 7 = YR 2 /n = 1480, from 
which we can extract the ratio of the Young's modulus Y 
to the bending rigidity n, given that the virus diameter 
is 2R = 60 nm |3j . The precursor capsid or "prohead" of 
HK97 is rather spherical, in contrast to the larger, more 
faceted mature infectious virus shown in Fig. ^] |2j| ■ It 
seems likely that this virus particle undergoes a buckling 
transition as it passes from the prohead to its mature 
infectious form. Indeed, the prohead shell is wrinkled 
or corrugated relative to the mature form It seems 
reasonable to regard this transformation as a change in 
the effective thickness d of the viral shell: d! of the pro- 
head goes to d < d' in the mature form. As discussed 
in Sec. Ill Al 7 = 12(1 — i> 2 )(R/d) 2 if we approximate the 
shell by a uniform isotropic elastic medium with Poisson 
ratio ^3 5] . Since 2R increases from 52 nm — > 60 nm as 
d decreases, it is plausible that 7 ~ {R/d) 2 rises and that 
the transformation from prohead — > head is accompanied 
by a buckling transition. 

Figure IT21 shows a similar fit to the yeast L-A virus, 
which yields 7 = YR 2 / k = 547. The diameter of the 
yeast virus is 2R = 43 nm, which leads to the conclusion 
that [F/K]yoast = 1-24 nm -2 . Note that |Y/«]hk97 = 
1.64 nm -2 , consistent with the arguments given in the 
Introduction that spherical viruses have roughly similar 
clastic constants. 

Our models for virus shells are based on two impor- 
tant assumptions. The first is the neglect of a sponta- 
neous curvature term j2(J in bending energies such as 
Eq. |5J Such a term might be significant if the viral build- 



ing blocks had a pronounced conical shape similar to, 
say, surfactant molecules in a micelle [26|, or labo- 
ratory cork stoppers. Several lines of evidence suggest 
that neglect of this term might be justified. Although 
certain virus scaffolding proteins (which can act as tem- 
plates for early phases of construction) do have a coni- 
cal shape, these are discarded in the mature icosahedral 
viruses of interest to us here [ij . It is hard to see a strong 
mechanism for precisely-defined hinge angles in very large 
viruses composed of many capsomers |29l \3(j . In vitro 
assembly experiments on the polyoma virus do produce 
spherical aggregates with 12, 24 and 72 pentameric units 
depending on conditions of pH, calcium concentration, 
etc [3l]], which could be accounted for by a spontaneous 
curvature term. Although we neglect spontaneous cur- 
vature here, it is certainly possible that the physics of 
small viruses (or the scaffolding proteins themselves) are 
influenced by such an energy, as has been explored by 
Bruinsma, Rudnick and Gelbart [32|. 

In general our work is more likely to be applicable to 
large viruses, for which simple continuum models can be 
justified. Icosahedral viruses are usually composed of a 
combination of 5- and 6-fold symmetric packing units, 
with the twelve 5-fold units centered on the vertices of 
an icosahedron (the polyoma virus SV40 Q with its 72 
identical pentamers is an exception). Because the strain 
energies which lead to bucklin g ex tend far from the discli- 
nations which produce them |13||. we would expect our 
results for large viruses to be insensitive to differences 
in the shapes of packing elements. Special packing el- 
ements at the twelve 5-fold sites could be incorporated 
into a disclination core energy. 

A second key assumption is our neglect of the osmotic 
pressure due to the confined DNA or RNA package of the 
virus Here we can appeal to an experiment. Earn- 
shaw and Harrison |3^| have compared the structure of 
phage lambda (P22) with its full complement of DNA 
to the structure of lambda mutants containing only 78% 
of the native DNA. Although changes in the details of 
DNA packing can be detected, the protein shell itself is 
unchanged. Thus, due either to DNA condensation or 
an exceptionally strong shell, the osmotic pressure of the 
DNA is insufficient to change the shape. Of course, the 
nucleic acid content of a virus could nevertheless play an 
important role in shell assembly [34|. 
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FIG. 11: Real space fit to a virus structure. The curve shows the root-mean-square deviation from the experimental virus 
shell of bacteriophage HK97 (full virus and cross section shown in the lower part of the figure) and the theoretically calculated 
shape. The best fit occurs at the minimum for 7 ~ 1480 and the corresponding shape is depicted in the inset. Shading indicates 
the distance of the shell from the center. 
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FIG. 12: Real space fit to Yeast virus L-A with inset and actual virus as described in Fig. 1111 The best fit occurs for 7 = 547. 
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